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Abstract 

The thermalization of a longitudinally expanding color glass condensate with Bjorken boost 
invariant geometry is investigated within microscopical parton cascade BAMPS. Our main focus 
lies on the detailed comparison of thermalization, observed in BAMPS with that suggested in 
the "Bottom-Up" scenario. We demonstrate that the tremendous production of soft gluons via 
gg — > ggg, which is shown in the "Bottom-Up" picture as the dominant process during the early 
preequilibration, will not occur in heavy ion collisions at RHIC and LHC energies, because the 
back reaction ggg — > gg hinders the absolute particle multiplication. Moreover, contrary to the 
"Bottom-Up" scenario, soft and hard gluons thermalize at the same time. The time scale of thermal 
equilibration in BAMPS calculations is of order ot~ 2 (In a s )~ 2 Q 'J 1 . After this time the gluon system 
exhibits nearly hydrodynamical behavior. The shear viscosity to entropy density ratio has a weak 
dependence on Q s and lies close to the lower bound of the AdS/CFT conjecture. 
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I. INTRODUCTION 
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The success of employing simple ideal hydrodynamics [l[ in describing the large values 
of the elliptic flow v 2 measured in Au+Au collisions at the Relativistic Heavy Ion Collider 
(RHIC) [2|, yj indicates that thermal equilibration of the produced quark gluon system occurs 
on a short time scale and that the equilibrium is maintained until the hadronization. It is 
of great interest to understand what mechanisms drive the system to equilibrium. While 
coherent quantum effects like color instabilities [4] may play a role at the very early stage 
when the system is super dense, perturbative QCD (pQCD) bremsstrahlung processes are 
essential for momentum isotropization of quark gluon matter [^[fl when the matter becomes 
more dilute due to the strong longitudinal expansion. 

The importance of pQCD bremsstrahlung was first raised in the "Bottom-Up" scenario 
which describes the thermal equilibration of a color glass condensate {J characterized 
by a saturation scale Q s . The main idea of the "Bottom- Up" scenario is that while the hard 
gluons with transverse momenta of order Q s degrade as the condensate evolves in space time, 
soft gluons with transverse momenta much smaller than Q s are populated due to pQCD 
99 ~^ 999 bremsstrahlung. Soft gluon production dominates the early stage of equilibration 
and a strong parametric enhancement of the soft gluon number has been predicted. Within 
a short time scale the soft gluon number becomes comparable to the initial number of 
hard gluons. As soon as the radiated soft gluons achieve thermal equilibrium and build up a 
thermal bath, the hard ones begin to loose their energy to the thermal bath and subsequently 
thermalize on a later time scale. A parametric time scale for overall thermalization is given 
by r th ~ a-^Q-' fl. 

Because color glass condensate (CGC) [8| is proposed as a possible initial state of the 
quark gluon matter produced in high energy heavy ion collisions, its thermalization is a 
highly interesting topic. For instance, the thermalization of an idealized form of color glass 
condensate was studied by Bjoraker and Venugopalan by [ll| solving the Landau trans- 
port equation. Serreau and Schiff also investigated the same topic in (l2| where they used 
the relaxation time approximation to simplify the collision term in the Boltzmann equation. 
A conclusion, which can be drawn from the two studies, is that pQCD gg — > gg collisions 
are not sufficient to achieve thermal equilibrium. Hence, the inelastic gg <-> ggg processes 
are needed as emphasized in the "Bottom-Up" scenario. 
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The role of the pQCD gg <-> ggg processes for thermal equilibration employing the color 
glass condensate initial condition is investigated in the present work for the first time within 
a full 3+1 dimensional transport calculation using the parton cascade BAMPS [kJ. The 
bremsstrahlung processes indeed lead to rapid thermalization. In this paper we focus on 
how thermalization occurs in BAMPS calculations. In particular, we investigate whether 
the "Bottom-Up" scenario is realized as the proper way to describe thermalization of CGC, 
which may be formed in heavy ion collisions at RHIC and LHC energies. Note that the back 
reactions of the bremsstrahlung, i.e., ggg — > gg processes, are consistently incorporated in 
BAMPS, whereas they are missing in the "Bottom-Up" scenario. This, which is shown 
in our results, leads to a different thermalization picture from the one suggested by the 
"Bottom-Up" scenario. 

This paper is organized as follows. In section [Til we present the parton cascade BAMPS 
and the setup. A boundary condition is introduced to mimic a one-dimensional (longitudi- 
nal) expansion with Bjorken boost invariance. In section II I II the detail on the CGC initial 
conditions is given. The importance of ggg — ► gg processes, which hinders the multiplication 
of soft gluons, is discussed in section HVl The numerical results based on BAMPS calcu- 
lations are shown in section |V] We compare the thermalization observed in BAMPS with 
the "Bottom-Up" scenario, determine the time scale of thermal equilibration for various 
saturation scales Q s and coupling constants a s and extract the shear viscosity to entropy 
density ratio, rj/s. A conclusion is given in section IVT1 

II. PARTON CASCADE BAMPS 

BAMPS is a microscopical transport model, which solves the Boltzmann equations for 
on-shell quarks and gluons using Monte Carlo techniques. In the present work only gluons 
are considered. The main feature of BAMPS is the implementation of 2 — > N and N — > 2 
processes in a consistent manner, which is based on the stochastic interpretation of interac- 
tion rates j^]. The interaction rates or the interaction probabilities are calculated locally in 
space where the phase space density of particles, f(p,x), is extracted numerically. BAMPS 
subdivides space into small cells, which are regarded as the local positions where interac- 
tions may occur. The smaller the cells the more local interactions can be realized. However, 
smaller cells contain fewer particles and thus lead to larger uncertainties in the extraction 
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of f(p,x). Therefore, we adopt the test particle method to amplify the (pseudo)particle 
density by a factor N test . The cross sections have to be reduced by the same factor to obtain 
the same physical mean free path 5]. 

In order to make comparisons with the "Bottom-Up" scenario we calculate the space 
time evolution of gluons in a tube with a radius of R = 5 fm. The transverse wall of the 
tube serves as a boundary to mimic one- dimensional (longitudinal) expansion. Gluons are 
simply reflected on the cylindrical wall. In the transverse plan a static "spider web" like 
cell structure is considered: The polar angle and the radial length squared r 2 are divided 
equally within [0, 2ir] and [0,i? 2 ], respectively. This division gives the same transverse area 
for all cells. For the numerical calculations we set A0 = 7r/4 and Ar 2 = 5 fm 2 . Longitu- 
dinally, space is divided in Az bins, which have the same width in the space time rapidity 
r) — j ln((i + z)/(t — z)). Ar] = 0.2 is set to be a constant for all Az bins. The initial gluons 
are put into rapidity interval [—3; 3]. Our setup is adequate to the assumption of Bjorken 
boost invariance 14], which is used in this study. 

Gluon interactions included in BAMPS are elastic pQCD gg —>■ gg scatterings as well as 
pQCD inspired bremsstrahlu ng qq «- > qqq. The differential cross sections and the effective 



matrix elements are given by 
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dq 2 ± (g 2 +mp 2 ' (1) 



|2 V s 2 I2g\\ 



M gg -* ggg \ = — j- 2 ■ o . 2 rjTTr ■ _ 21 9(fc ± A g -cosh|/) (2) 



2 (qi + m 2 ,) 2 k 2 ± [(k ± - q± ) 2 + 



in 



where g 2 = 4na s . q^ and k^ denote the perpendicular component of the momentum 
transfer and of the radiated gluon momentum in the center-of-mass frame of the collision, 
respectively, y is the momentum rapidity of the radiated gluon in the center-of-mass frame, 
and A g is the gluon mean free path, which is calculated self consistently [5J. A discussion of 
the present idealistic implementation of the LPM effect is given in jg]. 
The interactions of the massless gluons are screened by a Debye mass 

d 3 p 
(2tt)V 

where do = 16 is the gluon degeneracy factor for N c = 3. mo is calculated locally using the 
gluon density function f(p,x) obtained from BAMPS. The suppression of bremsstrahlung 
due to the Landau-Pomeranchuk-Migdal effect is taken into account within the Bethe-Heitler 
regime using a step function in Eq. (j2J). 



m 



D 



f d?p 

^ dGas JTw~v Nch (3) 



The interaction rates per particle are obtained [5j by 

R22 = n(v rc \a 2 2)2, (4) 
-R23 = n(v rc \a 2 3)2 , (5) 

2 { 8E1E2E3 ^3 ^ ^ 

for 00 -> 00, 00 -> 000, and 000 -> 00, respectively, where 
1 /" s / 4 da"^ 9g 

022 = -,\ dq 2 ± , (7) 



2! Jo ^ dqi 

am = hhl dT '^ dT '* l^™! 2 (2vr) 4 5 (4) (pi + - Pi - P ' 2 - p' 3 ) , (8) 
h2 = k\S dT ' ldT ' 2 \ M ^^ 2 ( 27r )^ (4) (Pi + P2 + Ps - p'i - P2) (9) 

where dT^ = d 3 p'J (2n) 3 2 E' i: \M ggg ^ gg \ 2 = \M gg ^ ggg \ 2 j 'da, s is the invariant mass for the 
interaction, v re \ = sj2E\Ei denotes the relative velocity of two incoming gluons and (-}2 and 
(•)3 symbolize the average within ensembles of incoming gluon pairs and triplets, respectively. 
For each gluon pair and triplet positioned in a cell unit with volume AV the transition 



probability within a time interval At is given 
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P22 = v vc ia 2 2-^, (10) 

At . . 

P23 = v^^—, (11) 

I32 At 

p ' 32 = mE 2 E3 av 2 ' (12) 

respectively, as derived directly from the transition rates (jlj"®- Note that R32 = 3i?23/2 for 
thermal equilibrium. The factor 3/2 indicates the ratio of the number of identical particles 
in the initial state of the 000 — > gg and 00 — > 000 interaction. 



III. CGC INITIAL CONDITIONS 

Jl 

For the initial condition a gluon distribution of a color glass condensate [8j is applied. 
The theory of a color glass condensate is given by the saturation picture, which assumes 
that the parton distribution in a hadron or nuclei saturates at high energies as a result of 
competition between QCD bremsstrahlung and annihilation processes. 

The CGC initial condition used in our simulations consists of gluons with Pt < Q s , which 
are produced by the nonperturbative part of the nucleus-nucleus interaction. Q s denotes 
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the saturation momentum, which is the typical momentum of g luons in the CGC. It is close 



181 ] . The color glass condensate 



to 2 GeV at RHIC and is expected to be 4 - 6 GeV at LHC 
is a state with high parton occupation number where th transverse momenta reach up to 
Q s , whereas the occupation number drops to for transverse momenta much larger than 
Q s . Initially, most gluons have transverse momenta close to Q s , whereas the longitudinal 
momentum of gluons in the central rapidity bin is approximately zero. 

For the initial gluon distribution of Color Glass Condensate we employ an idealized and 
boost-invariant form 



f(x,p) 



-s(y - vMQ 2 s - p 2 t) 



(13) 



a s N c r 

We take N c = 3 for SU(3). The factor c in ({TBI is the "parton liberation coefficient" which 
accounts for the transformation of virtual partons in the initial state into on-shell partons 
in the final state, as introduced in [20] . The value of c used in [llj was calculated for a 



SU(2) gauge 
of c ~ 0.4 m 



theory to be c = 1.3 [9j] 1 . New SU(3) gauge theory calculations yield a value 
lo| . which we employ for the following calculations 



The initial particle density in the CGC approach is given by lOj, LL3J 



1 dN 



N 2 



-Ql 



(14) 



irR 2 drj An 2 a s N c 

For the application of the Boltzmann equation, we need the phase space density to be smaller 
than unity. If phase space density is high, Bose enhancement factors should be considered 
in the collision integrals, which is not done in BAMPS model. 

The initial gluons are produced at eigentime r ~ ^ and the initial phase space density 
f(x,p) from Eq. (jT3l) is infinite due to the delta function 5(p z ) ~ Later the distribution 
in longitudinal momentum space broadens due to 2 — > 2 (or 2 — > 3) collisions and the 



occupation number becomes finite. Our cascade starts at time Tq 



a 3 N t 



■Ti where t,- 



At this time the parton distribution function in Eq. (|T3|) is still larger than unity. However, 
the same initial time has been applied in In the "Bottom- Up" picture at a time 

r ~ a~ 3 / 2 Qj 1 the distribution should become less than 1. We note that we do not employ 
the Bose enhancement factor (1 + /) within the Boltzmann collision terms. Hence, as long 
as / is larger than 1 we underestimate the rates. 



1 This value we used in our previous calculations 2l|, |22 1 . 
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In the following we present the results of simulations for Q s = 2, 3 and 4 GeV, i.e., 
energies relevant for RHIC and LHC. 



IV. SOFT GLUON PRODUCTION AND ANNIHILATION 

The basic assumption of the "Bottom-Up" thermalization picture is that soft gluons 
are continually emitted due to inelastic 2^3 bremsstrahlung processes, which increase 
in the soft the gluon number and leads to formation of a thermal bath. However, the 
annihilation processes, which are the back reactions of bremsstrahlung and missing are in 
the "Bottom-Up" scenario, will hinder the soft gluon production due to detailed balance. 
The annihilation processes become dominant, when the system is oversaturated. Using the 
CGC initial conditions in Eq. (fT4"l) we estimate whether a net production of gluons is possible 
during the subsequent space time evolution. 

Assuming a free streaming of CGC initial condition, the energy and particle densities are 
given by 



n t 



r 47r 2 a s N c s r 
Also assuming an instantaneous thermalization at time r t h with 



n ( T o)- = c— — -Q s -. (16) 
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e(r th ) = e th (r th ) = — T 4 (r th ) , (17) 
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one obtains 



The ratio 

w *h(n h ) _ 24 /J_\ ! A c -\ k T \ (19) 
n(r th ) _ UsJ s Qs th ( j 

estimates whether a net production or a net reduction of gluons will occur at the early 
times of the expansion. If the ratio ^th( r th)/ n ('H:h) is larger than 1, more particles should be 
produced. Fig. [T] depicts the ra th (r tll )/r7,(r tll ) ratio for fixed Q s = 3 GeV and various values 
of ClS db function of the thermalization time rth- All curves start at values smaller than 
1, which indicates that for the chosen parameters indeed the annihilation processes would 
dominate the early stage of equilibration. For an increase of the gluon number at early times, 
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FIG. 1: reth( r th)/ n ( T th) ratio in Eq. (fT9|) as a function of r t h for Q s = 3 GeV and a s = 0.1, 0.3, 
0.6, and 0.9, respectively. 

as predicted in the "Bottom-Up" scenario, the value of a s has to be much larger than 0.3, or 
Q s has to be chosen much larger than would be given at RHIC and LHC For small coupling 
the particle number will start growing if the time scale of thermalization is large. In the full 
calculation of the Boltzmann equation the gluon thermalization is more complicated than 
the simple consideration of a free streaming and a subsequent instantaneous equilibration. 
However, the behavior of the n t h(r)/n(r) ratio holds, as we will shortly see. 

The initial distribution of gluons is highly anisotropic in momentum space. Most gluons 
have a transverse momentum of order p t ~ Q s . Populations of the low (high) momentum 
gluons should be dominated by 2 — » 3 (3 — > 2) processes. Figure [2] shows the gluon p t 
spectra after one single 2 — >• 3 or 3 — > 2 interaction. While 3^2 processes increase gluon 
number in high momenta, 2^3 collisions lead to an enhancement of gluon number in soft 
momenta, which resembles the "Bottom-Up" scenario. 
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FIG. 2: Transverse momentum distribution after one single inelastic scattering. Initial condition 
is a CGC with a s = 0.3 and Q s = 3 GeV. 

V. RESULTS: THERM ALIZATION OF A CGC 
A. Gluon number 

Fig. [3] shows the gluon multiplicities per space time rapidity, dN/dr], at midrapidity 
(77 G [—0.1,0.1]) as function of time, which are obtained using BAMPS with the CGC 
initial conditions with a s = 0.3 and Q s = 2, 3, and 4 GeV, respectively. Their ratios to 
the initial gluon number are depicted in Fig. |H The way that thermalization proceeds 
within the parton cascade calculations does not resemble the way that has been advocated 
in the "Bottom- Up" scenario [7]. The strong parametric enhancement of the total gluon 
number at early times, as predicted by "Bottom-Up" scenario, is not observed in the cascade 
calculations. Instead, gluon annihilation occurs during the first 0.3 — 0.75 fm/c for Q s = 2 — 4 
GeV. This is clearly due to the 3 — ^ 2 annihilation processes and indicates that the initial 
CGC is oversaturated for the chosen values of a s and Q s . Figure [5] shows the interaction 
rate of both elastic and inelastic processes in the central space time rapidity bin. The rate 
of 3 — > 2 processes is initially significantly higher than that of 2 — > 3 processes, which leads 
to a decrease in the net gluon number at very early times. 

The gluon number begins to increase after t ~ 0.3 — 0.75 fm/c (depending on the value 
of Q s ) when the system is close to kinetic equilibrium and a quasi-hydrodynamical cooling 
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FIG. 3: dN/drj, at midrapidity (ij £ [—0.1,0.1]) as a function of time. Results are obtained using 
BAMPS for the initial CGC with a s = 0.3 and Q s = 2, 3, and 4 GeV, respectively. 
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FIG. 4: Ratio of the gluon number to the initial one in the central space time rapidity bin. 

sets in. Slow increase in gluon number at late times is consistent with the collision rates, 
i?23 ~ 2i?32/3, shown in Fig. [5l (In chemical equilibrium i? 2 3 = 2R 32 /3.) Assuming parton- 
hadron duality, the final gluon multiplicities, dN/drj ~ 700 — 1700 for Q s = 2 — 3 GeV, are 
equivalent to dN^/drj m 470 — 1100 for the total charged mesons, which are comparable 



with the RHIC data for Au+Au most central collisions at y/s = 130 - 200 AGeV [23j, [24]. 
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FIG. 5: Interaction rates in the central space time rapidity bin, obtained from BAMPS calculations 
for CGC with a s = 0.3 and Q s = 2, 3, and 4 GeV. 

As demonstrated in Fig. [2], 2 — > 3 collisions lead to an enhancement in the number of 
the soft gluons, whereas 3 — > 2 processes initially increase number of gluons with momenta 
higher than the saturation momentum Q s . We now study the changes in gluon momenta for 
various interaction channels included in the calculations. For this we define a soft momentum 
scale, which is p^ o{t < a s Q 2 s in the "Bottom-Up" scenario, to be p so f t = 1.5 GeV and a hard 
scale Phard = Qs- The medium gluons are denoted as gluons with p soft < p t < Phard- These 
definitions are in particular reasonable at early times since the longitudinal momenta are 
very small. They are different from the definitions in [7] , where the authors set all the initial 
gluons to be hard. Fig. O shows the net gluon production of each type in 2 — > 3, 3 — > 2, and 
2^2 processes as function of time. The results are in terms of the difference between the 
number of outgoing and incoming gluons of each type divided by the number of collisions. 
We see that 2 — > 3 collisions increase the soft gluon number with a loss of medium and hard 
gluons [see Fig. EJ^a)], whereas 3 — > 2 processes increase the hard gluon number with the 
loss of soft and medium gluons [see Fig. [0(b)]. Equivalently, 2 — > 3 processes transfer energy 
from the medium and hard to the soft sector, whereas 3^2 processes transfer energy from 
the soft and medium to the hard sector. Compared to 2 <-> 3 processes, 2^2 convert 
few medium gluons into soft and hard ones [see Fig. EJ^c)]. Common to all three collision 
types in Fig. [6] is a continuous degradation of the medium sector with the simplified CGC 
initial conditions. As the system evolves towards equilibrium, energy is transferred from the 
medium into both the soft and hard sector. 

The ratio of the numbers of the soft, medium and hard gluons to the total number of 
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FIG. 6: Net production of soft, medium, and hard gluons in various collisions. Results are obtained 
in the central space time rapidity bin from a calculation performed for a CGC with a s = 0.3 and 
Q s = 3 GeV. 

gluons is depicted in Fig. [7J The total gluon number is dominated by the medium sector 



s 

o 

a 

E 
-a 

<D 

E 
i= 
o 
to 



0.5 1 1.5 2 2.5 3 

t (fm/c) 



FIG. 7: Ratio of the numbers of the soft, medium and hard gluons to the total number. 

until 0.5 fm/c and then by the soft sector after ~ 1 fm/c. Contrary to the "Bottom- Up" 
picture, Fig. [7J shows that the soft gluon number increases over a long period of time at the 
cost of the primary "medium" gluons. To repeat, the production of soft gluons is effectively 
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hindered by 3 — > 2 processes and, thus, cannot exhibit a huge increase as predicted in the 
"Bottom-Up" scenario. 

B. Kinetic equilibration and momentum isotropization 

Fig. Ogives the transverse momentum spectra in the central space time rapidity and at 
various early times (up to 0.5 fm/c) obtained from BAMPS for initial CGC with a s = 0.3 
and Q s = 3 GeV. After the expansion starts, energy flows immediately into both the soft 
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FIG. 8: Transverse momentum spectra in the central space time rapidity and at various early 
times. The initial condition for the BAMPS calculation is a CGC with a s = 0.3 and Q s = 3 GeV. 

(pt < p S oft = 1-5 GeV) and hard momentum region (p t > Q s ) where the populations rapidly 
increase, as also seen in Fig. [7J Note that at 0.5 fm/c the number of soft gluons is of the 
same order as the number of harder gluons. 

From Fig. [H] we observe that the spectrum of high momentum gluons achieves an expo- 
nential shape on a short time scale and almost as quick as the soft gluons. However, they 
have different slopes. At t = 0.5 fm/c the entire spectrum is to a good extent in agreement 
with a thermal fit taking T ~ 0.67 GeV, which is indeed very close to the effective temper- 
ature of the system at this time [T e s(t = 0.5 fm/c) = E/3N = 0.6 GeV]. The transverse 
momentum spectrum achieves a thermal shape in hard and soft regions almost simultane- 
ously. The presence of a thermal bath of soft gluons seems not to be a necessary condition 
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for the equilibration of hard gluons. Again, this is different from the picture invoked in the 
"Bottom-Up" scenario. 

In 25] where the dynamics of SU(2) gauge fields in presence of an initial anisotropy 



in momentum space is studied, it is shown that the energy obtained from the particles 
by a Weibel-like plasma instability does not lead to an exponential buildup of transverse 
magnetic fields. Rather it is transferred into the ultraviolet modes via a rapid "avalanche" . 



This phenomenon, which was also discussed in [261 ]. is similar to what we have observed 
during the very early thermalization of CGC (see Fig. [8]). A more detailed study of such a 
highly nonlinear phenomenon is certainly needed. 

For CGC initial condition with a s = 0.1 and Q s = 3 GeV the rapid "avalanche" is again 
seen in from Fig. The number of hard gluons with an exponential shape increases on 
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FIG. 9: Same as FigJHJ The initial condition is CGC with a s = 0.1 and Q s = 3 GeV. 

the same time scale as in the case for a s = 0.3. This can be understood from the following 
consideration: Because at early times 7 32 [see Eq. fl9])] is roughly proportional to a 2 and the 
gluon density is inversely proportional to a s , the initial interaction rate, R32 ~ n 2 I 32 [see 
Eq. (jHJ)], is approximately independent of a s . 

Fig. HO] shows the transverse momentum spectra for a s = 0.3 at times larger than 0.5 
fm/c. The spectra are compared with thermal fits using temperature T as a parameter. The 
effective temperatures T = E/3N, which is extracted from the simulations at various times, 
are found to be indeed close to the values chosen for the fits. Thus, the transverse spectrum 
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FIG. 10: Same as FigJH]for later times. 

at t — 0.5 fm/c looks almost thermal. Moreover, the cooling of the system sets in, which 
makes the exponential spectra steeper at later times. This behavior is characteristic for a 
hydrodynamical expansion. 

Summarizing the results above, the thermalization of a CGC, studied in the parton 
cascade BAMPS, is characterized by the following facts to be opposed the "Bottom-Up" 
scenario: 

• Emission of hard (p t > Q s ) gluons due to 3 — > 2 processes clearly dominates the very 
early evolution. 

• No strong enhancement of total gluon number is observed. The total gluon number 
decreases slightly, until the system is nearly thermalized. 

• No thermal bath of soft gluons is built up at very early times. 

• Transverse momentum spectra achieves a thermal shape in the hard and soft regions 
almost simultaneously. 

C. Time scale of thermalization 

Next we extract the time scale when the system is more or less thermalized. A parametric 
dependence of the time scale, ~ a~ was given in the "Bottom-Up" scenario |7| . Fig. 
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[TT1 shows the effective temperature, T, and the scaled one, Tt 1//3 , as a function of time 
obtained with a s = 0.3 and Q s = 2, 3, and 4 GeV, respectively. The effective temperature 
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FIG. 11: (a) Effective temperature and (b) Scaled effective temperature. 

increases slightly at very early times due to gluon annihilation and then falls approximately 
with t -1 / 3 [see Fig. ITHfb)]. The same scaling is also found for the transverse energy per 
rapidity at late times, as seen in Fig. fl2T a). where dEx/dy ■ t 1 ' 3 is depicted. The behavior 
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FIG. 12: (a) Scaled transverse energy per rapidity and (b) scaled energy density, 
at late times corresponds to a one-dimensional ideal hydrodynamical expansion. Thus, we 
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determine the time scale of thermalization as the time, at which Tt 1//3 becomes a constant. 
The times extracted from Fig. [TlTb) read: 



UhiQs = 2 GeV) = 1.2 fm/c, t th (Q s = 3 GeV) = 0.75 fm/c, t th (Q s = 4 GeV) = 0.55 fm/c . 

(20) 

The value for Q s = 3 GeV is consistent with the time scale (0.5 fm/c) at which the transverse 
momentum spectrum becomes exponential (see Fig flQ|) . In addition, the values of t t h-Q s are 
almost equal for fixed a s and various Q s , which verifies the relation t th ~ Q^ 1 as predicted 
in the "Bottom-Up" scenario. 

Although the expansion at late times is described nearly by ideal hydrodynamics, the 
collision rates are not infinitely high (see Fig. [5]). Thus, the viscosity might be small 
but still nonzero (as shown in the next subsection). The effect of the nonzero viscosity is 
actually apparent in the scaled energy density, et 4 / 3 , depicted in Fig. H2T b). Although et 4//3 
is constant in time according to one-dimensional ideal hydrodynamics, et 4//3 still increases 
at later times in the simulations. 

For fixed Q s = 3 GeV and various a s the time scale of thermalization is extracted from 
Fig. [13] such that 

t th (a s = 0.1) = 1.75 fm/c, t th (a s = 0.2) = 1 fm/c, t th (a s = 0.3) = 0.75 fm/c, (21) 

and 

t th a s 13/5 = 0.0044 fm/c(a s = 0.1), 0.015 fm/c(a s = 0.2), 0.033 fm/c(a s = 0.3) . (22) 

The dependence of t th on a s proves to be considerable weaker in our cascade calculations 
compared to what is estimated in the "Bottom-Up" scenario. In Ref. jy] the authors found 
that the time scale of thermalization is the inverse of the total transport rate, which is 



proportional to a s 2 (lna s ) 2 T 1 28] (T ~ Q s in the present case). This scaling indeed 



holds for the thermalization times from our calculations: 

t th a 2 s (lna s ) 2 = 0.09 fm/c(a s = 0.1), 0.1 fm/c(a. = 0.2), 0.1 fm/c(a s = 0.3) . (23) 

Therefore, the time scale of thermalization is of order a~ 2 (lna s )~ 2 Q~ 1 , which is smaller 
than the "Bottom-Up" prediction. The quick thermalization observed here is consistent 
with the findings from the previous studies BE3.Q : The gluon bremsstrahlung favors 
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FIG. 13: Scaled effective temperature in simulations with fixed Q s = 3 GeV and a s = 0.1, 0.2, and 
0.3. 

large-angle radiation due to the LPM suppression, which is the reason for the dominance of 
the pQCD gg <-» ggg processes in thermal equilibration. 

Figure [H] shows the time evolution of the total gluon number per space time rapidity. The 



o 

II 




FIG. 14: Gluon number per space time rapidity in calculations with fixed Q s = 3 GeV and a s = 0.1, 
0.2, and 0.3. 

initial condition with a s = 0.2 and Q s = 3 GeV could be appropriate for Pb-Pb collisions 
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at the maximal energy at the LHC In this case the gluon number also decreases at the 
beginning. Thus, "Bottom-Up" thermalization might not be favored at LHC. 

Momentum isotropization, which is an important part of kinetic equilibration, can be 
described by the time evolution of (p 2 /E 2 ) [6j shown in Fig. [T5J where CGC initial conditions 
are used with a s = 0.3 and Q s = 2, 3, and 4, respectively. (•) denotes the average over gluons 




FIG. 15: Momentum isotropization. CGC initial conditions with a s = 0.3 and Q s = 2, 3, and 4, 
respectively. 

in the central space time rapidity. Due to expansion, (p 2 z /E 2 ), which is extracted in a finite 
spatial volume, saturates at a value that is slightly smaller than its equilibrium value 1/3. 
We fit (p 2 z /E 2 )(t) at time t using a relaxation ansatz 

n<)4 + (<|)W-i)ex P (-^) (24) 

where 9 xe y gives the time scale of momentum isotropization [6]. The fits are shown in Fig. 
H3 Using t = 0.5 fm/c, we find 6 rel = 0.85, 0.52, and 0.42 fm/c for Q s = 2, 3, and 
4 GeV, respectively. These time scales are smaller than thermalization time scales [see 
Eq. (l20l) ]. which indicates that momentum isotropization is completed before (nearly) full 
thermalization with quasi-ideal hydrodynamical expansion. 
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D. Ratio of the shear viscosity to the entropy density 



As already noticed in the previous subsection, e t 4 / 3 increases slightly at later times [see 
Fig. [T2Tb)]. which shows a deviation from ideal hydrodynamics. However, the agreements 
of T t 1//3 and dE^/dy ■ t 1 / 3 with ideal hydrodynamics [see Figs. [TTland [T2T a)] indicates that 
the shear viscosity (or better the ratio of the shear viscosity to the entropy density) is small. 
In the Navier-Stokes approximation, the diagonal elements of the stress tensor are given 



271 ] in the rest frame by 



( du- 1 -> \ 

T u = P - 2 V f ^ - -V • u\ - kV ■ u (25) 

where 7] denotes the shear viscosity, k the bulk viscosity, and P the pressure. We then obtain 

T +T —IT 

± xx~ ± yy - 1 - zz i nn\ 

7] = - — - , (26) 

2 (3 d z u z - V • u) 

S P — T — T — T 
k = x v (27) 

3 V ■ u 

For the system of massless gluons where e = 3P, the bulk viscosity vanishes. The flow veloc- 
ity u is expected to be approximately the same as given in . For an ideal hydrodynamical 
expansion u ~ (0,0, z/i) when the gluonic system is (nearly) thermalized. Thus, we obtain 

77 = l - (T xx + T yy - 2T ZZ ) (28) 

where T xx , T yy and T zz can be extracted from the numerical calculations. Note that the 
divergence of the flow velocity can be better extracted from the simulations because it 
relates to the transport collision rate, as derived in j^ . 

Due to large numerical uncertainties the entropy density it is difficult to extract s from 
the simulations. Therefore, the formula 

S = ~ l /eq0 °' X) /cq(P ' X) " 11 = ^ ~ U lll(A) (29) 

with feq(p, x) = Xdce~^ is used, which is applied to systems in kinetic equilibrium. Here n 
is the gluon density and A = n/nth denotes the gluon fugacity where n t h = dcT^/ir 2 is the 
gluon density in thermal equilibrium. The entropy density calculated in Eq. ( 1291) is, thus, 
larger than the true value in the simulations. 

Figure [TBI shows the r]/s ratio in the calculations with a s = 0.3 and Q s = 2, 3, and 4, 
respectively. Before 0.5 fm/c the values of rj/s are not reliable because at early times the 
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FIG. 16: Ratio of the shear viscosity to the entropy density. CGC initial conditions with a s = 0.3 
and Q s = 2, 3, and 4, respectively. 

gluon system is far from equilibrium and, thus, Eqs. fl28l) and fl29l) are not valid. From 0.5 
fm/c on the t]/s ratio is nearly constant and has a weak dependence on Q s , rj/s ~ 0.15, 
which is exactly the same as that obtained in full 3+1 dimensional BAMPS calculations 
with a s = 0.3 and minijets type initial conditions for Au+Au collisions at RHIC energies 



291 ] . This verifies that the rj/s ratio determines the behavior of the late dynamics and, thus, 



depends only on the coupling a s , but not on initial conditions. Moreover, the pi s ratio 



obtained is small and close to the lower bound from the AdS/CFT conjecture [30j. The 
smallness of the rj/s ratio corresponds to the efficiency of the pQCD gg <-> ggg processes in 
thermal equilibration, because the rj/s ratio is inversely proportional to the total transport 
collision rate and the transport collision rate of gg <-» ggg processes is 6 — 7 times larger 



than that of gg — ► gg collisions 28] 



VI. CONCLUSION 



Using the parton cascade BAMPS, we have studied the thermalization of potential color 
glass condensates, which might be appropriate for the initial conditions of high energy heavy 
ion collisions. The main emphasis is put on the comparison of the thermalization observed 
in our calculations to that in the "Bottom-Up" scenario. We found that several aspects 
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of the real thermalization might be different compared to the "Bottom-Up" scenario. The 
difference arises because the back reactions of bremsstrahlung, 3^2 processes, play a 
significant role. They are completely absent in the "Bottom-Up" scenario. 

First, the radiation of gluons is hindered by 3 — > 2 processes according to detailed balance. 
Therefore, for the idealistic form chosen for CGC, the total gluon number will increase unless 
a s Q s is very large, which is not realistic at RHIC and LHC We showed that for realistic 
initial conditions the total gluon number decreases early in the expansion. An exorbitant 
increase in soft gluons, as predicted in the "Bottom-Up" scenario, is not possible. Thus, no 
thermal bath of soft gluons will be built up. 

Second, thermal equilibration of soft and hard gluons occurs roughly on the same time 
scale due to 2 — ► 3 and 3^2 processes, respectively. The energy flows into both the 
soft and hard sectors at the same time, which is potentially similar to the phenomenon of 
"avalanche" as observed in the field isotropization driven by the plasma instability. This 
behavior contradicts the "Bottom-Up" picture where soft gluons form a thermal bath and 
thermalize first whereas hard gluons loose energy to the thermal bath and, thus, thermalize 
later. 

Finally, the time scale of thermalization is determined for various values of a s and Q s . It 
spreads from 0.55 fm/c to 1.75 fm/c for Q s = 2 — 4 GeV and a s = 0.1 — 0.3. In agreement 
with the "Bottom-Up" scenario, the thermalization time proves to be proportional to Qj 1 , 
however, its proportionality to a~ 13 ^ 5 is not seen, but is much weaker: r t h ~ (a s In as) -2 ^ 1 . 

After being thermalized the gluon system shows quasi-hydrodynamical behavior: The 
cooling due to expansion is observed in the steepening of the transverse momentum spectra. 
To see how viscous the system is, we extracted the ratio of the shear viscosity to the entropy 
density and obtained r)/s ~ 0.15 for a s = 0.3. The rj/s ratio has a weak dependence on Q s 
and is close to the lower bound from the AdS/CFT conjecture. Thus, the considered gluon 
system acts as being strongly coupled. 

The quick thermalization and the smallness of the rj/s ratio observed in the present 
calculations with the CGC initial conditions are consistent with the findings from the pre- 



vious studies 



28 



29| using the Glauber-type minijets initial conditions. This 



demonstrates that independent of the chosen initial conditions, the pQCD bremsstrahlung 
processes (and the back reactions) dominate the dynamical equilibration and then keep the 
system behaving like a nearly perfect fluid. The higher order processes such like ggg — > ggg 
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and gg «-> gggg will certainly lead to a larger total transition rate, however, their contribu- 



311 ]. Further investigations are needed to quantify 



tions are suppressed by higher order of a s 
these contributions. 

The CGC initial conditions presented in this paper are idealistic. More realistic initial 
conditions in high energy nucleus-nucleus collisions can be obtained by considering high mo- 



mentum jet like partons(minijets) [32j]. This will modify the thermalization picture presented 
in this paper: ggg —>■ gg interaction may become less dominant at the early stage of the 
thermal equilibration, however, it would not affect our conclusion that thermalization in the 
hard and soft sectors proceeds on the same time scale. A study will be presented elsewhere. 
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